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Abstract 

In engineering applications of fuzzy logic, the main goal is not to simulate the way the experts really 
think, but to come up with a good engineering solution that would (ideally) be better than the expert's 
control, [n such applications, it makes per feet sense to restrict ourselves to simplified approxi m ate 
expressions for membership functions. If we need to perform arithmetic operations with the resulting 
fuzzy numbers, then we can use simple and fast algorithms that are known for operations with simple 
membership functions. 

In other applications, especially the ones that are related to humanities, simulating experts is one of 
the main goals. In such applications, we must use membership functions that capture every nuance of 
the expert’s opinion; these functions are therefore complicated, and fuzzy arithmetic operations with the 
corresponding fuzzy numbers become a compu tational problem. 

In this paper, we design a new algorithm for performing such operations. This algorithm is applicable 
in the case when negative logarithms — log(/i(r))of membership Lunctions /r(x]are convex, and reduces 
computation time from 0(n 2 )\.o 0(n log(n)) (wheren is the number of points xat which wc know the 
membership functions /i(xj). 


1 Formulation of the Problem 

Depending on the goal, applications of fuzzy logic can be naturally divided into two classes: 

• Engineering applications like fuzzy control in which fuzzy logic is used as a tool for achieving a certain 
goal: a better (smoother and safer) control of a car, a better heating, etc. In such applications, 
the expert’s knowledge described by fuzzy rules is used not to simulate the way experts solve these 
problems, but to design better control strategies, 

• Applications to humanities (psychology, linguistics, etc. ) in which fuzzy logic is used to describe and 
simulate the human behavior, the human decision-making processes, etc., and thus predict the way 
humans will react in dilTerent situations. 

In both types of applications, we have to deal with fuzzy numbers r, i.e., quantities whose values we do not 
know precisely, and instead, we only have expert (fuzzy) knowledge about these values. This knowledge is 
usually described in terms of membership functions /j r (x) that assign to every real number x the expert’s 
degree of belief p r (x)£ [(), 1] that the actual (unknown) value of the quantity r is equal tox. 

The formalism (membership functions) is the same, but, depending on the application, wc treat these 
membership functions differently: 

• In engineering applications, we do not need to describe the exact opinion of the experts, because wo 

are going to improve this description (b-j some fine-tuning) anyway. Therefore, it is quite sufficient to 
use membership functions that approximately describe export’s opinions. To simplify computations, 
usually, the simplest approximations are used, most often triangular or trapezoid membership functions 
(see, e.g., [2]). 
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• In humanities applications, if we use oversimplified approximations to membership functions, we will 
end up having very crude models of human behavior For such applications, we, therefore, need accurate 
descriptions of membership functions, and these descriptions can be very complicated. 

1.1 Fuzzy Data Processing and Fuzzy Arithmetic Operations: if We Must Use 
Precise Membership Functions, We Have a Computational Problem 

Fuzzy data processing. We want to use the expert (fuzzy) knowledge about the values f'l , . . . . r u of 
some quantities to predict the value of some quantity r that is related to r;. I n this paper, we will consider 
the simplest case when “related "means that we know the exact form of the dependency r= /(r r„) 

between and r, and the only uncertainty in r is caused by the uncertainty in the values of r, . 

For example, when we formalize the expert’s opinion about possible candidates for a position, we may 
know that this opinion depends on the values of it characteristics r i of the candidate, we have expert (fuzzy) 
knowledge about the values of r >, and wo know that the final opinion depends on the total evaluation 
r = wi ■ r\ + . + w n r n with known weights Wi. 

In such situations, we must transform the fuzzy knowledge about the valuesr,- into a fuzzy knowledge 
about r = f(r j, . ,,r„). This transformation is called fuzzy data processing. 

Fuzzy arithmetic operations. In the computers, usua]ly : only elementary arithmetic operations { + i 
■, /) are hardware supported. Therefore, every data processing algorith m written in a high-level program- 
ming language is parsed, i.e., represented as a sequence of elementary arithmetic operations. For example, 
computing an expression Xj {*2 + 13) is decomposed into two steps: Computing £2 + X3 and multiplying 
the result by xi- 

In view of this decomposition, in order to implement, an arbitrary data processing algorithm with fuzzy 
inputs, it is sufficient to be able to apply elementary arithmetic operations o = +, — , to fuzzy numbers. 

The formulas for these operations come from the extension principle (see, e.g., [3]): In particular, if we use 
an algebraic product a - b as a fuzzy analogue of fz. we arrive at the following formula for t = r o s: 

/i,(x) - sup (p r {y) (1) 

y,z yoz—x 

[n particular, for o = +, we have 

p t {x) = s»p(/i r (.y) ■ p,(x - y )). (2) 

v 


For simple membership functions, fuzzy arithmetic operations are computationally easy. For 

example, if we use Gaussian membership functions 

Pr(x) — exp{(x - a r ) 2 /{a r f), 
p,{x) =exp((x- a, f !((?., f), 

then (2) leads to a Gaussian membership function for t: (x) — exp ( (x - a t ) 2 ) with 

a r (cr r )~ 2 +a a ((T s )-' 2 

at M " 2 + ( a ,)~ 2 

and ((r ( ) -2 = (oy) -2 + (r t , )-2 [3 ; 5 ] 7 hese are computationally verv simple formulas to implement. 

There are simple formulas for several other cases (see, e.g., [3] and references therein). 

For complicated membership functions, fuzzy arithmetic operations are computationally com- 
plicated. When wc cannot use approximating simple expressions, then we cannot use simplified formulas 
that stem from the use of these expressions, and therefore, wc have to use the formula (2), This formula is 
straightforward, so, we can simply use it to compute /^(x). To find out how long it would take to compute 
Pt(x), let us estimate tile number of computational steps that arc required to compute pt(x). 

Of course, in reality, wc can only know the values of p r (x) and p>{x) for finitely many values x. Let, us 
denote the total number of such values byn. In this case, it is reasonable to compute only n values of pt(x). 
For each of these n values, according to the formula (2), we must find the largest of n products. Computing 
each product takes 1 elementary computational step, computing the largest of n numbers requires that we 
don — 1 comparisons. So, the total number of computation steps that needs to be done to compute one 
value of p. t (x) is 2 n - 1 = 0(n). 
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If we have n parallel processors at our disposal, then we can use each processor Lo compute its own value 
of p t (x) and thus, compute all these values in linear time. 

In many real-life situations, however, we only have one computer. In such situations, to compute all n 
values of the desired membership function pt(x), we need 0(n 2 ) computational steps. 

The more accurately we wish to represent the expert’s opinion, the larger n we need to take. For large 
n,0(n 2 ) is too long. Can we perform fuzzy arithmetic operations faster? 

In [4], an approximate algorithm is given that performs arithmetic operations with fuzzy numbers in time 
0(n log(n)). 

1.2 What We Are Planning to Do 

In this paper, we design a new fast algorithm that computes the precise value of the resulting membership 
functions in 0(nlog(rt)) time. 

This algorithm is applicable when the negative logarithms - log(^t(a;) ) of the membership functions 
p(x) are convex. This class of membership functions includes many important classes such as Gaussian 
membership functions. 

2 Fast Addition of Fuzzy Numbers 

2.1 Main Idea 

Let us describe, step-by-step, how we can simplify the problem of computing the sum of two fuzzy numbers. 

First simplification: reformulation in discrete terms. We only know the membership functions p,(x) 
and p.,(x) in finitely many points, and usually, these points are of the type ar>= i Ax. In this case, the 
formula (2) takes the following form: 

ti = ma x(rjS,_j), (3) 

where we denoted ti = p t (i Ax),r,~ p r (i . Ax), and Si — p,(i Ax). 

Further simplification: reducing multiplication to addition. The formula (3) can be simplified even 
further if wc recall that the equality t= r • s is equivalent to ‘T = R. + S, where T = - ln(<), R = - ln(r), 
and S = - in(s). In view of this equivalence, and taking into consideration the fact that - in(z) is a strictly 
decreasing function, we can reformulate the formula (3) as follows: 

Ti =miri (llj + Si-j ), (4) 

i 

where we denoted 7\ = - ln(f,-),fJi = - ln(rj), and S, = - ln(s;). We will describe how, given the two 
sequences Ri and S , , we will be able to compute the elements T; fast. Then, if we know the values r,- = 

p r (i-Ax) and s,- = p s (iAx), we will be able to compute the values Ri and Si, compute T, = - ln(<,), and 
then reconstruct the desired values ti = pt(i Ax) as ti = exp(— Ti). 

How to compute the formula (4)? 

Final simplification: a local criterion for the maximum. Fora given i, when does the sum £T = 

R.j + Si-j attains its minimum’? If it does attain the minimum for some j, this means that the value of this 
sum for this particular j is not larger than the values of this sum for j - 1 and for j + 1: -i anc * 

J2j 5; 52j+i • ^ we denote Dj = then these two inequalities take the form 

Dj < 0; D i+ , > 0. (5) 


We can use binary search to find the desired j. Since the function - ln(^ r (x)) is convex, the sequence 
Rj is also convex, and therefore, the differences R.j- Rj_ y are monotonically increasing with j. Similarly, 
the differences )- are strictly decreasing with j. Therefore, the difference Dj = ( Rj- Rj-y) - 
(Si-j- Si-^j- 1 )) is increasing with j. 

Hence, we can find the desired value j that satisfies the condition (5) by using binary search: This will 
be an iterative process on which, on each step, wc will have lower and upper bounds for the desired value j. 
We start with the lower and upper bounds that encompass all possible values of j. Then, on each iteration, 
we: 
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• take a midpoint m = (lower + upper) div 2 between the current lower and the upper bounds; 

• compute D m for this midpoint, m, and 

• compare the resulting value D m with 0. 

Depending on the result of this comparison, wc do the following: 

• If D m = 0, then, due to the monotonicity of the sequence D m , we have D m + 1 > D m = Of ■ - e -f 
D m+ i > 0. Hence, this m satisfies the condition (5). Using monotonicity of Dj , one can easily show 
that in this case, 

- either m is the only value for which (5) in true (in which case, it is the only possible minimum of 

2Zm)’ 

" or Z? ; =0not only for j= m, but also for several values of j that are neighboring tom, in which 
case, there arc several minima with exactly the same value of ■ 

In both cases, the value of for the midpoint m is the desired minimum. 

• If D m > 0, this means, due to nionotonicity of the sequence Dj , that j < m. In this case, we can take 
m as the new value of the variable upper. 

• Similarly, if D m < 0, this means, due to monotonicity of the sequence Dj , that m < j. In this case, 
we can take m as the new value of the variable lower. 

This algorithm takes 0(n log(n)) steps. On each iteration of the binary search, we reduce the size in 
half. In k iteration, we go down from n to <nj 2 k possible values. When nj 2*<I, we are down to a single 
point, and thus, wc have localized the desired j. The inequality n/2 k < 1 is achieved when k ss log 2 (n), so, 
we need 0(log(n)) points to find the desired j and thus, to compute the desired value of 7; for this particular 
i. 

To compute the values of Ti for 71 different i’s, we thus need n 0(log(n))= 0(n iog(n» computational 
steps. 

2.2 Resulting Algorithm 

GIVEN: the values p,(r) and jt,(x) for n equally spaced values*: = r’ . Ax. 

ALGORITHM: 

• First, for each of n values x,-, we compute the values Ri = - ln({i r (x t )) and 5',- = - (*:'))• 

• For each i, we: 

_ apply binary search to find the index j for which the non-decreasing sequence 

Dj = ( Rj- Rj - 1 ) - (Hi-j- Si -7 j ) passes from the non-positive to non-negative values; 
compute T, as/?j+ Gi-j for this very j ; 

‘ compute Rt(xi) asexp(— T,). 

3 Algorithms for Other Arithmetic Operations 

3,1 Subtraction 

To compute t — r -s, we can represent it as t= 7*+ ( _s )- Since we know the membership function /.i s (x) for 
s, we can easily compute the membership function /<_ s (x) for -s as /i_ s (z)= /i,(— x). Then, wc can apply 
the above algorithm to compute the desired membership function for l = r - s = r + (— s). 
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3.2 Multiplication 

If Lhe quantities r and s both take only positive values, then, to compute r s, we can use the formula 
r s = exp(ln( r) + ln(.s)) : 

• From the membership functions for r and s, wo can easily compute the membership functions for ln(r) 
and ln(s) as W n ( r )(x) = p r (ln(z)) arid /i| n i»)( a; ) = Mln(x)). 

• Applying the algorithm presented above, we compute the membership function /<i n (t) for ln(f) = In(r)+ 
ln(s), 

• Finally, from we compute ^t(.v) as [i t {y) = p.| n ( t )(exp(y)). 

3.3 Division 

Division t= r/s can he expressed as t = r ( 1/s). So, to divide two fuzzy numbers, we can use the following 
algorithm: 

• First, we compute the membership function for 1/sas \iyj s [x.)- 

• Then, wc use the algorithm for multiplication to compute the membership function for 

t = r ■ (l/.s) = r/s. 


3.4 Computational Complexity 

For all these operations, the major part is computing the sum of fuzzy numbers that takes 0(n log(n)) steps. 
Therefore, the computational complexity of computing the difference, product, or ratio of two fuzzy numbers 
is also 0(nlog(n)). 

4 What If A t-Norm (^-Operation) Is Different From Algebraic 
Product? 

4.1 Fuzzy Arithmetic Operations: Case of a General t-Norm 

For an arbitrary & — operation /&(rt, b), the extension principle for addition leads to the following formula: 

l*t(x) = sup fsAnr{y), H, ( X - ?/)). (3) 

4.2 Strictly Archimedean t-Norms and Reduction to the Case of Algebraic 
Product 

Idea, [t is known (see, e.g., (3]), that if an &— operation satisfies some reasonable conditions, then it can 
represented i n the form 

fa{a,b) = t/'-’CVKa) ■ ip(b)) (4? 

for some strictly increasing function ip : [0, 1] —+ [0, 1] (&— operations that satisfy these "reasonable" condi- 
tions are called strictly Archimedean). 

Since the function ip is strictly increasing, the value /&(/ir (y), /i,(x - xj)) is the largest ifT the value 
t/))) is the largest, so. 

Hl*t(x)) = sup ^’(U.(^r{y),H,{xy)). (5) 

Y 

From (4), we conclude that xp(f^(fir(y),li,(x — y)))=xl>(fi r {y))-ii , {l-‘A x ~ y)). Therefore, (7) can be rewritten 
as: 

'/■'(MaOl'sup iifb{y)} {x - y )). (6) 

y 

If we denotery(a:)= \p(n r (x)),v s {x)= 0(p r (i)),andi/ t (.s)= ip(fi t (x)), then this formula will take the form 

i/ t (x) = sup(i/ r ( y)v»( x ~ Y)) ! (7) 

Y 
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which is exactly like the formula (2) that we already know how to compute fast. From Vt(x)= tl’(pt(x)), we 
can compute fit(x) by applying an inverse function Pt{x) = (ft(x)). 

So, to compute p t (x),wc can apply the following algorithm: 

Algorithm. 

• For every x , compute u r (x) = lf(p r (x))axidl/ s (x) = lp(fl s (x)). This takes O(n) steps. 

• Apply the algorithm (described in the previous section) to i/ r (x) and u„ (x); this algorithm will take 
0(n log(n)) computational steps and return nu t (x). 

• Apply the inverse function to i/ t ( x), resulting in ij>~ 1 (v t (x)) . This is done value-by-value, 

so, for O(n) values of x, it takes 0(n) steps. 

Computational Complexity. The resulting algorithm requires 

o(n) + 0(nlog(n)) + 0(n)= O(nlog(rt)) 


computational steps. 

4.3 Other Arithmetic Operations 

For other arithmetic operations with fuzzy numbers (-, . . /), we have a similar reduction to the case of 
algebraic product that leads Lo similar 0{n log(n)) algorithms. 
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